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TECHNICAL NOTE 2857 

A THEORETICAL METHOD OF ANALYZING THE EFFECTS 
OF YAW-DAMPER DYNAMICS ON THE STABILITY OF AN AIRCRAFT 
EQUIPPED WITH A SECOND-ORDER YAW DAMPER 
By Albert A. Scby and Ordway B. Gates, Jr. 


SUMMARY 


A method is described for investigating the effects of the dynamic 
response of an autopilot on the stability of an aircraft-autopilot com- 
bination. The method is based on a study of the constant-damping curves 
obtained in a plane defined by varying tvro of the autopilot parameters. 
The dynamics of the autopilot are assumed to be describable by a second- 
order differential equation. The effects on the system stability of 
varying the gain, natural frequency, and damping ratio of the automatic 
damper are investigated, since these parameters determine the dynamic 
response of the automatic damper. 

The method is applied to the analysis of the lateral motion of an 
airplane equipped with a second-order automatic yaw-rate damper. For 
any condition of the aiiylane, em optimum combination of values of auto- 
pilot natural frequency and damping ratio are shewn to exist for any 
given gain or required damping. A simple, analytical expression is 
derived for obtaining a close approximation to these optimum points' by 
ignoring the effects of the aperiodic characteristic modes of the air- 
plane. The assumption that these aperiodic modes may be neglected in 
considering the effect of the yaw damper on the Dutch roll oscillation 
is used in all the subsequent analysis. Expressions are derived for 
the maximum damping obtainable under various conditions. For any given 
natural frequency and damping ratio of the autopilot, excessive auto- 
pilot gain will always cause the autopilot oscillatory mode to become 
unstable . 

Finally, the problem of designing an efficient yaw damper which 
will improve the damping of the Dutch roll oscillation for vaxious 
flight conditions of an airplane is considered. A simple method of 
design is lUtistrated by applying it to three fll^t conditions of an 
airplane . 

Calculated motions, based on the assimiption of three degrees of 
freedom for the lateral airplane motion, are presented. They agree 
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with the resialts obtained from the constant- damping- curve analysis when 
the aperiodic airplane modes are neglected. 


INTRODUCTION 


Recently, a great deal of interest has been shown in the use of 
automatic stabilization devices for Improving the damping of the lateral 
(Dutch roll) oscillation of aircraft designed to travel at transonic and 
supersonic speeds. The analyses of such automatic stabilization systems 
may be divided into two classes. In one kind of analysis the effects of 
various types of autopilots on aircraft stability are considered, eind 
usually the autopilot is asstnned to be an idealized system with no lags. 
In the other kind of analysis the effects of the dynamic response of a 
particular type of stabilization system on the stability of the aircraft- 
autopilot combination are considered. 

Some investigations of the effects of various types of idealized 
autopilots on aiiplane lateral stability are reported in references 1 
to 3* References 1 and 2 are theoretical analyses, whereas in refer- 
ence 3 experimental results are compared with theoretically predicted 
effects . 

Analyses of the second kind, in which the effects of certain types 
of dynamic response in a given stabilization system are considered, have 
involved various approaches to the problem of determining the stability 
of the complete system. The we3J.-known frequency- response analysis has 
been used in many papers (for example, in ref. 4). Applications of 
Nyquist’s stability criterion, which was originally developed for feed- 
back amplifiers, have also been used. For example, in reference 5 
Nyquist’s criterion is extended to systems with constant time lag in 
the feedback circuit. 

When the conditions which will instire a given amount of stability, 
rather than j\xst neutral stability, axe sought, the use of constant- 
damping curves has been found convenient'. The method is described in 
reference 6, and an example of its application is furnished by refer- 
ence 7* In reference 8 a semigraphical method is developed for obtaining 
the conditions which insure neutral stability of an airplane-autopilot 
system when a constant time lag in the autopilot is assumed; also, a 
procedure is indicated for using the method of constant-danqjlng curves 
tQ determine the conditions which insure a given amount of damping. 

This method is described in detail in reference 9* 

In the present paper the constant-damping- curve analysis is applied 
to an airplane-autopilot system in which the autopilot dynamics are 
represented by a second-order differential, equation. Generally, the 
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actiial frequency response of most autopilots (especially yaw dan5)ers) 
is tetter approximated wten a second-order differential eqiiation is 
used to describe the dynamic characteristics of the autopilot than wlien 
the constant-time-lag response is assumed. Physically, this assvmqption 
ing)lie6 that the autopilot may be represented as a damped oscillatory 
system. The present paper examines in some detail the effects on the 
aiip>laiie“a’Utopllot stability' of variations in the gain, - natural fre- 
quency, and damping ratio of -the autopilot. • • ■ 


SYMBOLS 


A, B 


Cl = 


coefficients of second-order differential equation for auto- 
pilot dynamics (see eq. (l)) 

wing span, ft 




Cl 

Cn 

Cy 


Sc- 




Sc- 


2V 

Sc? 

Cl = 

‘r > rb 


Kx" 

trim lift coefficient, W cos 7/qS 
rolling-moment coefficient. Rolling moment/ qSb 
yawing-moment coefficient. Yawing moment/ qSb 
lateral-force coefficient. Lateral force/qS 

■ ST 

SCy 

% = ^ 

^ . 

SCy 

% . 




5Cn 

S5 


^Ur - 


_ ^n 


P Sp 


c - 


Cn_ — 


SCn 


D 


time-derivative operator, d/dt 
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F(D) 

Fo(D) 

fn(R) = 

i = ^ 

K 

Ko 

% 

Kz 

%Z 

m 

n 

P, Q 

Pq^ Qo 

p = 

q 

R 

r = Dt 
S 

^1/2 

t 


characteristic polynomial of aiip)lane-autopilot system 
characteristic polynomial of airplane alone 

: d.°F(R) 
dR° 


gearing ratio of second-order autopilot, deg^deg/sec 

gearing ratio required for an ideal no-lag rate damper to 
provide a given amoxint of damping, deg^deg/sec 

nondimenslonal radius of gyration in roll about longitudinal 
stability axis 


nondimensionaJL radius of gyration in yaw about normal 
stability axis 

nondimenslonal product-of-ineortia parameter 


mass of airplane, slugs 
an integer 

coefficients of a ■ quadratic factor of characteristic equation 
of airplane -autopilot system 

coefficients of Dutch roll quadratic (see eq. (lO)) 


dynamic pressure, Ib/sq ft 

real part of characteristic root, sec"^ 

wing area, sq ft 

time for amplitude of oscillation to danrp to one-half its 
original value, sec 

3 

time, sec 
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V 

W 

3 

y 

5 

p 

03 


“A 


03o 


steady-state velocity, ft/sec 
wei^t of airplane. It 

angle of sideslip, radians tinless otherwise specified 
fli^t-path angle, radians 

deflection of control stirface, deg or radians 
damping ratio of second-order autopilot 

I 

relative-density factor, m/ pSh 
air density, slugs/cu ft 

angle of hank, radians unless otherwise specified 

angle of yaw, radians unless otherwise specified 

angular frequency (always referred to simply as "frequency") ^ 
radians/sec 

frequency of Dutch roll oscillation, radians/sec 
natural frequency of second-order autopilot, radians/sec 


[]]]^(D) transfer function of aiiplane 

p(D) transfer function of autopilot 


Subscripts: 
cr critical 

max maximum 


ANALYSIS 

Preliminary Discussion 


The dynamics of an autopilot used for aircraft stabilization can 
often be represented by a differential equation of the type 

D^6 + ADS + B5 = KBX 


( 1 ) 
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Here 5 is the deflection of the control surface actuated by the auto- 
pilot, and X is the con^onent of airplane motion to which the auto- 
pilot is sensitive. In general, X may he any combination of the air- 
plane degrees of freedom or their time derivatives, or both. When B 
is positive, as it must be for stable autopilots, let A = 2^0^ and 

B = tOo^. Then equation (l) becomes 

D% + 2^o5qD5 + = Kc£>q^ (2) 

Although A and B are often more convenient to use for purposes of 
analysis, the equivalent parameters ^ and cUq are more meaningftil 

physically. Most of the results of this paper will therefore be pre- 
sented in terms of these parameters. The parameter £Dq is the natural 

(undamped) frequency of the autopilot system and the parameter ^ is 
the damping ratio, that is, the ratio of the actual damping of the 
system to the critical damping. The parameter K which appears in 
equations (l) and (2) is var-iously caJJLed the amplification factor, the 
gain, or the geajring ratio. An autopilot which may be represented by 
equation (l) or (2) is often called a second-order autopilot, since the 
dynamics of the autopilot are represented by an expression involving 
time derivatives of the control deflection up to and ’ including second- 
order derivatives . 

The coefficients of the characteristic equation of an airplane 
equipped with a second-order autopilot are functions of the stability 
derivatives emd mass characteristics of the airplane and of the three 
autopilot parameters K, cSq, and ^ (or K, A, and B) . These coef- 
ficients are often called the stability coefficients. If the stability 
derivatives and mass characteristics of the airplane are known for a 
given fli^t condition, and if one of the autopilot parameters is 
assigned some reasonable value, then the stability coefficients of the 
airplane- autopilot system are functions of the remaining two autopilot 
parameters only. CvLrves of constant dang)ing of the characteristic modes 
of the total system may therefore be drawn in the plane defined by con- 
sidering these two parameters as independent variables. Methods of 
obtaining such cirrves, with paJrtlc^ll£Lr application to oscillatory modes, 
are discussed in references 6 and J. From an examination of the constant- 
damping curves the ranges of values of the two independent autopilot 
parameters which will provide a given amount of damping to the oscilla- 
tory modes of the airplane -autopilo'E^ system can be determined. Moreover, 
a detailed study of these constant-damping curves yields a considerable 
amount of insist into some of the fundamental properties of the motions 
of aiip lanes equipped with second-order automatic stabilization systems. 
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Application to Second-Order Yaw-Eate Damper 

. Eqiiations of motion .- The method has been applied to the analysis 
of the effect of a yaw-rate damper on the lateral stability of an air- 
plane. When the dimensional time-derivative operator is used, the 
lateral equations of motion are 

2PbD3 - ^ CypP + 2p^,Dt - I = I Cy (3) 


- ^ + 2pi^Vilr - I I 


I I 


£ C 5 


V" 


(^) 


“ ^ - I I Ci^Dt + - | ^ ^ (5) 


In these equations the assumptions are made that Cj^ = Cy^ = 7 = 0 

and that the only result of the control deflection 5 is a yawing 
moment. The equation of motion of the yaw-rate damper, written in 
the form of equation (2), is 

D^8 •+ 25o)qD6 + 00q% = (6) 


The con5)lete characteristic equation for the eqioations of motion 
(3) fo (6) is a sixth-degree equation. In fact, if Fq(D) be the usual 
fourth-order characteristic polynomial, for the lateral degrees of free- 
dom described by eqimtions (3) to (5) (see ref. 2), then the complete 
characteristic equation for the airplane-autopilot combination is 


F(D) = (d2 + 2^0)oD + cDo^)fo(D) - ^ - I CY^)(2p-bKxV - 


~ J C7 d) - ^ CjPi 

2 b ^ y -h3 ^ n 


= 0 


(7) 
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The cliaxacteristic modes obtained from this equation are usiially two 
oscillations and two aperiodic modes. In most cases the two aperiodic 
modes and one of the oscillations may be associated with the lateral 
motion of the airplane, whereas the second oscillation may be associ- 
ated with the motion of the autopilot system. 

Choice of variable autopilot parameters for constant-damping 
curves .- If all the parameters in equations (3) to (^) are known for 
an airplane in a given fli^t condition, the effects of K, and ^ 

qn the lateral oscillatory stability of the combined system can be 
determined by varying these parameters in equation (l) , as is shown in 
appendix A. One of these parameters must be fixed if the stability 
bomdaries are to be plotted in a plane as described in appendix A, aind 
therefore the relative convenience of fixing each of the three param- 
eters should be considered. Since the dynamic characteristics of the 
autopilot may be expressed in terms of its frequency response, a study 
of the effects pf the three autopilot parameters on this frequency 
response should provide some insist into the question of which two 
parameters should be varied simultaneously. 

The transfer function of the stabilizing autopilot may be obtained 
from equation ( 6 ) : 


[5/D3p(D) = 


d2 '+ 


( 8 ) 


Therefore, the expression for the autopilot frequency response is 


[5/Di0p(hD) 


K 


1 - 


CDp 


CD 

<»o 




The freqxiency response of the autopilot is obtained from equa- 
tion ( 9 ) by plotting the amplitude Rp(cD) and phase angle 0p((D) of 

the complex nimiber [ 5 /D^p(ia}) against cd. The phase angle is inde- 
pendent of K, whereas the magnitude at any value of o) is proportional 
to K. Thus the phase-angle curve and the shape of the magnitude curve 
both depend only on ^ and co^^, and the gearing ratio simply acts as an 
amplification factor on the magnitude curve. From this point of view, 
it would seem desirable to select reasonable values of K and allow ^ 
and cDq to be the variable parameters in equation (T)* A study of the 
stability boundaries in the ^fJ^j-plane woiild then show the effect of 
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veirying ttie shape of the autopilot frequency response on the stability 
of the total system, whereas the variation of K would show the effect 
of varying the amplification factor. The stability boundaries were 
therefore first calculated in the ^ojQ-plane, and some of the effects 
of varying K on these ciirves were investigated. 

Since the gearing ratio K is in many respects the most iurportant 
parameter determining the stability of the total system, stability bouncf- 
arles were also obtained in the Kaig-plane for fixed values of ^ . This 
method of plotting the stability boundaries has two distinct adveintages. 
First, the two most impo2rtant parameters, K and cUq, are allowed to 
vary. Second, the fixed parameter ^ is known to be between 0 and 1 
for stable oscillatory autopilot systems; therefore, the effect of 
varying ^ on the stability boundaries in the Kh>o-plane can be deter- 
mined fairly easily by choosing several values of ^ which will span 
this range. 

Constant-damping curves with gearing ratio fixed; three degrees of 
. freedom .- In order to investigate the effect of the second-order yaw 
damper on the stability of the system, values of the airplane parameters 
were inserted into equations (3) to (5) to correspond to a cruising 
flight condition for the airplane described in table I. The Dutch roll 
oscillation for this flight condition has a period of 1.30 seconds and 
T1/2 of 2.60 seconds. 

In order td investigate the type of constant-damping curve which 
appears in the ^ojo-plane (or AB-plane) for constant K, a value of K 
, was chosen which would give good dan5)ing if used in a perfect propor- 
tional ya.w danger that has no inertia or damping. The value K = O.O86 
degree of rudder deflection per degree per second of yawing velocity 
was chosen, which wotild make Dutch roll oscillation equal 

to 0.75 second. The constant-dao^jing curves in the AB-plane were then 
drawn for this value of K by using the equations given in appendix A, 
and are shown in figures 1 to 5- The frequencies of the modes on a 
given cuive vary from zero to infinity, and typical values are shown on 
the individual ciirves. 

As pointed out in appendix A, it is simpler to obtain the curves 
first in terms of A and B and then substitute the more significant 
parameters ^ and cSq* However, for a preliminary investigation of 
the general types of curves and the manner in which they change as R 
and CO (the damping and frequency parameters) are varied, continued 
use of the AB-plane is more convenient because negative values of B 
may be considered, whereas the parameters ^ and cOq have an obvious 
physical significance only for positive values of B. For this reason, 
the typical cxirves of figures 1 to 5 are drawn in the AB-plane. These 
figures eire presented and discussed primarily to familiarize the reader 
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with the general types of curves which occur, so that the later com- 
parison of the curves in the ^oiQ-plane will not he too confusing. The 
practically significant portions of the curves are essentially the same 
in either plane. 

Figure 1 presents the zero-damping curve, which is continuous and 
crosses itself to fom a loop. The loop has no practical significance, 
since it merely defines the region -of autopilot parameters for which 
there are two unstable oscillations. The hatch marks indicate the 
boundary of the region in which all modes are stable. In all cases the 
B = 0 eixis is the boundary at which an aperiodic mode becomes tinstable. 

The curve of greatest practical importance is that defined by the 
damping of the airplane without the autopilot; that is, the curve for 
T 3_^2 = 2.6o seconds. This curve, shown in figure 2, is discontinuous. 

When CD approaches the aiiplane frequency (the frequency of the Dutch 
roll oscl3J.ation of the aiiplane alone) , both A and B become infi- 
nite in magnitude. These infinities are caused by the vanishing of a 
factor in the denominator of the expressions for A and B. At the 
airplane frequency a new branch of the constant-damping curve is started. 
The region in. the AB-plane in which there are no oscillatory modes that 
have less damping than the aiiplane al.one (indicated by hatch marks) is 
bounded by this new branch of the curve. The dashed line shows the 
value B = Qq = 23.84, where Qq Is the constant coefficient of the 

Dutch roll quadratic, which is written as D^ + PqD + Qq. The signifi- 
cance of this value will be discussed subsequently. 

Figure 3 shows a typical curve for a damping somewhat greater than 
that of the aiiplane alone. It has two points of discontinuity and 
three branches. Thus, the airplane-damping curve shown in figure 2 is 
a critical curve, separating the continuous curves for less damping 
from the doubly discontinuous curves for greater damping. The two 
critical values of <d are very close, so that the branch of the curve 
in the negative-B region corresponds to a very small range of a>. 

This branch of the curve is of academic interest only. The part of the 
curve in figure 3 which forms the significant boundary is again the 
curve in the upper right quadrant of the plane, defining a region of 
better damping indicated by hatch marks. 

The curve in figure 4 is for T^/2 =0*75 second, the damping 

which the airplane-autopilot oscillation would have if the autopilot 
had no lags. This curve has another branch at large negative values 
of B and negative A, but this branch is not shown because, as has 
been mentioned, it is of academic importance only. . (it should be 
noted that the significant portions of all these discontinuous curves 
start at values of cd near the airplane frequency.) The curve for 
T 2/2 = 0*75 second also represents a critical damping veilue, since 
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for larger values of damping the curves again become continuous and 
form loops, as shown in figure 5 for = O.60 second. In the 

following discussion the damping corresponding to = 0*75 second 

and the associated constant -damping curve will be called the critical 
damping eind critical curve. 


These curves will be discussed in greater detail. For the present 
it will suffice to note that, for any gearing ratio, greater damping 
can be provided by a second-order autopilot than could be provided by 
a perfect proportional autopilot with the same gearing ratio, provided 
that the proper values of ^ and o>o are chosen from a closed- loop 
region of the type shown in figure 

Equivalent-oscillator concept .- Appendix A shows that, in deter- 
mining the effect of the yaw damper on the stability of the Dutch roll 
oscillation, the airplane is represented by an equivalent oscillator 
whose period and damping are those of the Dutch roll oscillation. 
Explicitly,' the assimiption is made that the three-degree-of -freedom 
equations of motion of the airplane, given by equations (3) to (5)^ 
may be replaced by the single equation of motion 





+ P^Dt + 



C 6 - C 

^2 V -b2 


( 10 ) 


The physical interpretation of this assumption is that the effect of 
the aperiodic modes on the Dutch roll stability is small. The quadratic 
• P 

eqiiation D + PqD + Qq = 0 yields the complex characteristic root • 
corresponding to the Dutch roll oscillatory mode. 

The analysis is greatly simplified by this equivalent-oscillator 
concept, and it will be seen that this simplified analysis gives adequate 
results. Except where otherwise specified, the equivalent oscillator is 
assumed to represent the airplane in all the subsequent discussions. The 
curves in figiires 1 to obtained from the three-degree-of -freedom 
analysis, were drawn primarily as a check of the accuracy of the 
equivalent-oscillator approximation . 

In figure 6 the constant-danping c^lrves calculated by the equivalent- 
oscillator analysis for representative values of danping are presented in 
the ^cuQ-plane. In interpreting these curves it is essential to remember 
that two oscillations axe present for most of the significant points in 
the ^cOq - plane. One of these oscillations may generally be associated 

with the airplane and the other with the autopilot. Therefore, at every 
point in the most significant regions of the ^tOQ-plane two danping curves 
must cross each other. 
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All the curves previously ohtaiiied from the three-degree-of- freedom 
analysis are presented in figure 6, and in addition the boundary of 
eqiial roots and the curve for Tp/2 = 0*38 second are shown. Each of 
the constant-damping curves starts on the boundary of equal roots at 
(D = Q. Before the geneiral discussion of figure 6 is presented, the 
significance of the curve for = 0.38 second will be explained. 

Maximum damping with fixed gearing ratio .- As mentioned pi^eviously, 
the copstant-damping curves for hi^er damping than the critical damping 
are continuous and form loops defining the regions which insure greater 
damping. Since the loops become smaller for larger values of damping, 
it is reasonable to assume that for some value of damping the loop will 
become vanishingly small - only a cusp in the curve (see ref. 9 ) • This 
c\irve then corresponds to the maximum damping which can be obtained for 
the Dutch roll oscillation with a second-order yaw-rate damper for the 
given gearing ratio. Thus the cusp point may be considered the optimum 
point in the -plane for a given gearing ratio, and this point deter- 
mines the shape of the autopilot frequency response which will give the 
hipest damping to the Dutch roll oscillation for a given amplification 
factor K. 

These optimum points may be obtained by a rather simple algebraic 
analysis because they correspond to a double oscillatory root of the 
characteristic equation, which is a queirtic equation when the equivalent- 
oscillator concept is used. The derivation of the optimum-point char- 
acteristics is presented in appendix B. The damping and frequency of 
the double oscillatory mode are ass\nned to be given by the characteristic 
eqmtion 


(d^ + FD + q)^ = 0 


( 11 ) 


As shown in appendix B, the value of Q may be obtained by solving the 
quadratic equation, for a given gearing ratio. 


(" " = 0 ( 12 ) 

For positive values of K the larger real root in equation (l2) is used. 
With this value for Q, P may be obtained from the expression 

^ (Pq -I- ClK)Cp - PqQo^ 

2«o(e - 9o) 


(13) 
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These values of P and Q, when used in equation (ll)^ give the charac- 
teristic root (and therefore thq period and damping) corresponding to the 
hest-damped Dutch roll motion obtainable with the given gearing ratio. 

In order to find the autopilot psirameters A and B (or ^ and cdq) 
which yield this maximum damping, these values of P and Q are sub- 
stituted into the following expression’s for A and B: 


A = 2P ^ Pq 


( 1 ^) 


Qo 


(15) 


The maxi mum damping for the .airplane under discussion for 
K = 0.086 deg^deg/sec was obtained in this manner and was found to 
.correspond to Tjl/ 2 = 0*38 second. As can be seen from figure 6 , this 
curve does have a cusp at the optimum point in the -plane. 


Since the high-damping curves in the three-degree-of -freedom 
analysis also form loops, the same type of analysis , can be used to 
obtain the cusp point in the ^coQ-plane. The calculation of the maximum 
damping would be much more complicated, however, fqr the three-degree- 
of- freedom case. To show that the maximum damping as calculated by the 
equivalent-oscillator analysis -is an adeqtiate approximation to the maxi- 
mum damping for the conplete airplane, the curve for T^/2 = 0.38 second 

was drawn for the three-degree-of -freedom case also. 


The comparison of the curves for T^^2 =0.38 second for the two 

cases is shown in figure ?• lu this figure the significant portions of 
other typical damping ciirves are silso shown. For practical purposes, 
the equivalent-oscilJLator analysis is an adequate approximation to the 
three-degree-of -freedom analysis in determining the required constant- 
damping curves, including the maximum-damping point. It mi^t also be 
noted that the critical damping is not exactly the same in the two cases 
(^ 1/2 “ 0.75 second for the three-degree-of-freedom case and 

Tj^ 2 = 0*73 second for the eq'ui valent one-degree-of -freedom case^. 


Discussion and interpretation of damping boimdaries in -Plane . - 

Figure 7 shows that there is little change in the stability boundaries 
when the lateral motion of the airplane is represented by the equivalent 
oscillator. This result means that the real characteristic modes of the 
latereil motion can be neglected in calculating the effect of the auto- 
matic yaw damper on the Dutch roll oscillation. The only fundamental 
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difference occurs in the airplane-damping curve, ’^ 1/2 = 2.6o seconds. 

For the equivalent oscillator, this ctirve is discontinuous at a value of 
o>o equal to the aiiplane frequency, 4.8 radians per second. The slope 
of the curve becomes Infinite at this value of cDq, and the portion of 
the ciirve hounding the higher- damping region begins here at ^= 00 . in 
the three-degree-of -freedom analysis, the curve for T 1/2 = 2.6o seconds 
has a slope which becomes very large at cDq = 4.8 radians per second 
but remains finite. Thus the difference in the curves is negligible 
for practical piarposes. Since the significant portion of the airplane- 
damping curve begins at csq = 4.8 radians per second for any positive 
gearing ratio, the damping of the equivalent oscillator cannot be 
improved with a positively geared yaw-rate damper which has a natural 
frequency less than the frequency of the oscillator itself. Figure 7 
shows that for practical pirrposes the same statement can be made for 
the actual airplane. 

Certain general effects of varying the autopilot damping and natural 
frequency on the stability of the system may be observed in figures 6 
and 7* The area of interest is the rou^ly rectangular region indicated 
by the hatch marks in figure 2 in the AB-plane, in which both oscillatory 
modes have more damping than the airplane alone. The value of B indi- 
cated by the dashed line in figure 2 corresponds to the airplane fre- 
quency (cuo = 4.8 radians per second) . If ^ is fixed at any positive 
value and is increased, figure 6 indicates that the system damping 

increases to a maximum value at some value of cDq and then drops off, 
approaching the critical damping as csq approaches infinity. Similarly, 
at any fixed natviral frequency greater than the Dutch roll frequency, if 
the value of ^ is increased from zero the damping reaches a maximum at 
some value of ^ and then drops off, approaching the airplane damping 
as ^ approaches infinity. Thus, for fixed ^ there is an optimm cDq, 
and for fixed (Oq an optimum The best of all these points is the 

maximum-damping point, which is obtained by the simple calculation pre- 
viously described. 

Importance of oscillation frequencies in interpreting constant- 
damping curves .- From stability considerations alone, the regions defined 
by the hatch-marked portions of the curves shown in figures 1 to 5 deter- 
mine the values of autopilot ' parameters which guarantee at least the 
Indicated amount of damping. However, points may be chosen outside a 
given region which still seem to give an airplane motion that is as well- 
damped as that for points in the region. For exan^le, figure 8 shows 
the motions obtained with autopilots having the characteristics defined 
by points 1, 2, eind 3 in flg\ire 6. All motions shown in this paper were 
obtained from a Reeves Electronic Analog Computer, by use of the three- 
degree-of-freedom equations of motion. The fact that these motions check 
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the predictions made from the eguivalent-oscillator stability boundaries 
confirms the adequacy of the equivalent-oscillator analysis. 

Points 2 and 3 on the good portion (loop) of the curve for 
T2/2 =0.60 second, but point 1 is on the intersection of the curves 

for = 0-^ second and T2/2 = 2.6o seconds and is outside the 

good region defined by the loop. However, the actual airplane motions, 
represented by the sideslip and roll motions in figure 8, are very 
similar for all three autopilots, and the lightly damped mode is impor- 
tant only in the rudder motion for the case corresponding to point 1 

(fig. 8 (a)). Actually, the effect of the li^tly damped mode can be 

seen in the sideslip motion of figure 8 (a), but it is almost negligible. 
Although point 1 in figure 6 corresponds to = 2.60 seconds, as far 

as the ai2plane motion is concerned this autopilot would seem to give as 
good damping as autopilots whose characteristics fall in the loop of 
= 0.60 second. 

In order to understand why the lightly damped mode corresponding 
to point 1 in figure 6 has practically no effect on the airplane motion, 
the frequencies of the modes which are predicted at this point must be 
considered. Since for the most significant points in the ^o>Q-plane there 
must be two characteristic oscillations, each of these points is actually 
a crossing point of two constant-damping curves, as can be seen for 
point 1 in figure ,6. The general trends of the frequencies along the 
damping curves are shown in figures 1 to 5» Along the final portion of 
each curve (that is, the portion which approaches the ^ = 0 axis), the 
frequencies correspond to the autopilot frequency. At point 1 in fig- 
ure 6 the mode on the curve for T1/2 = 2'.60 seconds has a frequency 
03 « 10 radians per second, a vaJ.\ie which corresponds to the autopilot 
frequency (as can be seen from fig. 2), whereas the better -damped mode 
is the airplane mode with 03 5 radians per second (see fig. 5)* 

Since, the frequency of the autopilot mode is approximately twice that 
of the airplane mode, the effect of the corresponding lightly damped 
rudder oscillation on the airplane motion is sma.ll, for the airplane 
cannot follow such rapid oscillations. 

Consideration of the frequencies which occur at points along the 
damping ciirves is thus seen to be important in attempting to predict 
the type of motion which woiild be obtained with autopilots whose char- 
acteristics Eire deteimlned by these points. This frequency effect is 
brou^t out even more strongly by figure which shows the motions 
corresponding to points U and 5 in figure 6. In this case both points 
are on the zero-damping curve. At point 5 both modes have approximately 
the same frequency. However, at point h the autopilot frequency is the 
neutrally 'damped one, with 03 w 13*5 radians per second, and the airplane 
mode lies on one of the well- damped loop curves, which crosses the ciirve 
for point 4. The motions shown give the results predicted 
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by the analysis. As shown in figure 9 (a), corresponding to point k, the 
neutrally dan^jed hi^-frequency inidder motion has little effect on the 
well-damped airplane motion. On the other hand, for point 5 the effect 
of the neutrally damped mode, which is close to the natural airplane 
frequency, is the dominant neutrally damped airplane oscillation shown 
in figure 9 (b) . Autopilots which cause very poorly damped control 
motions, however, would be vinsatisfactory from a practical point of 
view even if their effect on the airplane damping were satisfactory. 

Effect of varying K on curves in ^o>o-plane.- The effect of geeiring 
ratio"'may be obtained by considering the effect of varying K on the 
curves in the ^oDq - plane. By using the equivalent-oscillator analysis, 

the critical dan5)lng is found to be Rcr = “^(^o Since the air- 

plane damping is -^o, the critical damping becomes the airplane damping 

as K vanishes. This result is, of course, necessary, since zero 
gearing ratio implies no autopilot. As K is made smaller, the whole 
set of loop curves in figure 6 tends to move to tho left, since the 
critical c\arve approaches the airplane damping curve. Conversely, as K 
increases, these curves move to the right. Also, the loop corresponding 
to any given damping larger than the critical damping must expand as K 
increases, since the given damping comes closer to the cfitical damping. 
Therefore the given loop approaches the infinite loop asymptotic to the 
critical-damping curve. Physically, this result simply means that as 
the gearing ratio is increased there is a larger range of values of ^ 
and (Uq for which a given damping lairger than the critical damping may 
be obtained. 


A clearer idea of the way in which the set of loop curves moves in 
the ^(Do-plane as K varies is obtained by investigating the vgiriation 
of the msLXimtmi-damplng point as K varies. The position of the maximum- 
damping point is itself important, because it is the optimum combination 
of ^ and cOq ^”0^ veilue of K; but, also, since this point is a 
kernel which is surrounded by all the loops, the motion of this point 
gives a clearer idea of the motion of any loop as -K varies. The 
desired variation may be easily obtained by inserting values of K into 
equations ( 12 ) and (l 3 ), and using the resulting values of P and Q 
in equations (l 4 ) and ( 15 )* However, a simpler and clearer method is 
shown to be possible in appendix B, wherein a valiie of damping is assumed 
and solutions are found for the values of K, A, B, and Q which will 
make this the maximum damping. This procedure clearly gives the smallest 
magnitude of K with which the desired damping may be obtained with the 
second-order automatic yaw damper, and the associated values of A and 
B then may be considered as "optimum" values for the given airplane and 
desired damping. In this method equation (13) is replaced by 


P = 


-2R = 


1.386 

Ti/2 


r 


(16) 
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where R and 
then gives A 


T 1/2 correspond to the desired damping. Equation (l4) 
immediately, and Q may he obtained from the exp>ression 


« - «o * - Po) 


( 17 ) 


Equation ( 17 ) will’ give two values of Q, one greater them Qq 
and the other less than Qq. By using these values in equation (15), 
the corresponding values of B are obtained, and from equation (13): 


K = 


2PQ - PqB - AQq 


CiB 


( 18 ) 


The smaller value of Q- obtained from equation ( 17 ) gives a value 
of B corresponding to a Value of- cDo lower than the airplane fre- 
quency and results in a negative K. The larger value of Q gives a 
value of (Dq higher than the airplane frequency and positive gearing. 
Figures 10 and 11 show the variation of the optimum points. Figure 10 
shows the cirrve on which the points lie in the ^o>o-plane. Figure 11 
shows the gearing necessary to obtain any given damping as the maximum 
damping (when the autopilot characteristics are the optimum ones for 
that gearing) . 

The point corresponding to K = 0 on each curve is significant 
only as a limiting point, since zero gearing implies no ■ autopilot . 

Since the double oscillatory mode corresponding to the maximum damping 
must approach the Dutch roll mode as K approaches zero, the optimum 
point in the ^flOQ-plane approaches the values corresponding to the Dutch 
roll mode. Increased damping of the system can be obtained with either 
positive or negative gearing. The positive-K branch of the locus of 
optimimi points given in figure 10 lies in the range of values of cOo 
higher than the airplane frequency. Thus, as mentioned previously, 
second-order yaw dampers with positive gearing must have a value of coq 
higher than the airplane frequency in order to improve the damping. 

Also, the variation of the optimum points on the positive-K branch of 
figure 10 shows that, if increased damping from the autopilot is sought 
by increasing the gearing ratio, the natural frequency and damping ratio 
of the autopilot should generally be Increased simultaneously. This 
fact can be of considerable practical importance, as will be brou^t 
out more clearly in subsequent discussion. 

Comparison of advantages of positive and negative geeiring .- If 
negative gearings are used, o)q values lower than the aiip)lane frequency 
must be used in order to improve the damping. The possibility of using 
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a negatively (reverse) geared second-order rate autopilot to improve the 
damping seems rather surprising, since for a perfect proportional rate 
autopilot, reverse geaxing would simply decrease the effective 0^^, and 

therefore decrease the damping. The use of reverse gearing is made 
possible by the phase relations introduced between the airplane and 
rudder motions by the dynamics of the second-order autopilot. Actually, 
this possibility is no more surprising than the fact that this type of 
autopilot will Improve the danping with positive gearing only when the 
autopilot natural frequency is greater than the airplane frequency. 
Clearly, this restriction arises' from the same type of phase-relation 
requirement. 

The use of a negatively geared yaw damper would seem advantageous 
because of its properties in a steady turn. For constant yawing velocity 
the negatively geared yaw damper deflects the rudder in a direction to 
maintain the turn, whereas the positively geared yaw damper must be 
overridden, either by the pilot or by the boost system. However, certain 
objections to the use of negative gearing in an automatic damper for use 
with an airplane actually make such use impractical. 

The main objection to the use of negative gearing is based on the 
fact that the equivalent oscillator represents the airplane only in a 
given fli^t condition. At different flight conditions the character- 
istic airplane oscillation has different values of damping and frequency; 
therefore, the airplane is represented by a different equivalent oscil- 
lator at each flight condition. To design an automatic damper for one 
fli^t condition only is impractical, since this automatic damper may 
have a harmful effect on the damping at some other fli^t condition. 
Autopilot characteristics must therefore be obtained by some compromise 
method which will improve any practical fli^t condition. Now, it can 
be shown that the regions of improved damping in the ^cnQ-plane for nega- 
tive values of K are loops resembling a reflection in the A = 0 axis 
of the unstable loop shown in figure 1. These loops must lie in a rela- 
tively narrow range of cDq values, since they are confined to values 
of (Do lower than the airplane frequency. Moreover, for a given magni- 
tude of gearing ratio, the loop for any damping is much smaller for nega- 
tive K than for positive K. Figure 11 shows that the loops for nega- 
tive values of K break down at much smaller values of damping than the 
loops for positive values of K of the same magnitude. If these small 
loops in the ^CD^-plane are drawn for a desired amount of damping for two 
extreme fli^t conditions with different natural frequencies, the possi- 
bility of their intersecting in a region of the ^(Do-plane which would 
give the desired damping to both flight conditions is relatively small. 

For positive geeirlng, on the other hand, improved damping can be 
obtained up to the critical damping for an infinite range of cDq starting 
at a value of cDq somewhat greater them the natural airplane frequency 
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in the particular fli^t condition. There is, therefore, an infinite 
range of cDq values which will insure at least the critical dan^jing 
for any number of fli^t conditions. The minimum valize of 03 q neces- 
saiy is somewhat greater than the highest natural frequency of any of 
the possible flight conditions of the airplane. 


Maximum damping for any gearing ratio .- The hipest maximum damping 
and the corresponding value of K can be obtained by determining where 
the modes corresponding to the optimum points in figures 10 and 11 become 

nono sc dilatory. This condition will occur for = 4Q. By using this 
qondition in equation ( 17 )^ where the positive sign corresponds to posi- 
tive K and the negative sign to negative K, the following expressions 
are obtained: 

W ^ + 1 j (K > 0) (19) 


■max 


= 2 « 


2 + 




(K < 0) 


( 20 ) 


Since Pg is small for li^tly damped airplanes, the positive gearing 
ratio gives a hi^er -value. For example, for Pg = 0, equation (19) 
gives a value approximately six times as large as equation (20) . 
Actually, these limiting values for the damping are of only academic 
interest as far as applica.tion to the airplane is concerned, since they 
are so large as to be far above any required damping. 


Constant-damping curves in KLOg -plane with damping ratio fixed.- In 

order to obtain a more complete picture of the effect of varying gearing 
ratio on the stability of the system, constant-damping curves were 
obtained in the King-plane with ^ fixed; As shown in appendix A, it 
is necessaory to solve a quadratic equation for the cDq -values, which 
may then be substituted into an expression for K, as follows: 


(Og - - flD^)cDg^ - 2 ^cDq(Pq + 2R)(o)^ + R^) +03^- 


(Og + 2 PgR + 2 R^)((d^ + R^) = 0 


(21a) 
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K = J(2^C0o + Pq + - 

L 


[jPo + 2 R)cDo2 + 2^oJo(Qo + 

(21-b) 

For the first calciolatlons, ^ =0-3 was asstuned as a reasonahle 
value of daB 5 )ing ratio for the autopilot. A set of curves was obtained 
for this value of ^ . The effect of varying . ^ on the curves in the 
KiOo -plane was then investigated hy obtaining several typical curves for 
^ = 0.6 and ^ = 0.9- 

The zero-damping curve for ^ =0.3 (shown in fig. 12) indicates 
that better-than-neutral damping can be obtained for any value of osq 
and for positive or negative gearing. !Fhis is true for any damping up 
to the airplane damping- The boundaries for less than the airplane 
damping are of no practical interest, however, and figure 12 is pre- 
sented only for completeness. 

Figure 13 shows the airplane damping boundary, which is a simple, 
continuous curve. This curve alone does not indicate clearly the. region 
that defines points which give better damping than that for the air- 
plane without yaw danger. However, the axis K = 0 must be part of the 
boundary also, since K = 0 Implies no autopilot, which means that the 
airplane has its original damping. In order to verify that the region 
defined by the hatching in figure 13 1b the good region, a curve was 
drawn for a sll^tly greater dan^jing (^ 1/2 ~ 2-50 seconds). This curve 
is shown in figure l4 and confirms the fact that the region insuring 
damping greater than that of the aiiplane is as shown in figure 13 . 

V 

The results of figxare 13 confirm several of the previous conclu- 
sions concerning negative gearing which were obtained from figures 10 
and 11. A relatively small region is present in the Kro-plane in which 
improved dan 5 )lng can be obtained with negative K, and this region is 
confined to frequencies less than the airplane frequency. 

Figure l4 shows that the regions of negative K which will give 
improved damping are loops in the Kuo-plane. Thus, the maximum damping 
for negative K may be obtained from the cusp point corresponding to 
the breakdown of these loops for any values of ^ . Since ^ is con- 
stant, a partic\ilarly sin^ple expression can be obtained for R,n„y. 

If ^ and cDq are used instead of A and B, equations (l4) and (l5) 
may be used in equation ( 17 )^ and the result is 


PqR + 3R^) + r(2Qo + 3PoR + 
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03o = ■ ^ ^ (K > 0) (22) 

= ~z (K < 0) (23) 

1 + t 


for the values at the cusp points. These values are given hy the 
two curves in figure 10. The maxlm\na damping for negative K at any 
value of ^ may then he obtained hy using equation (23) in equa- 
tion (l4) to yield 

If the values of Pq and Qq for the Dutch roll of the airplane under 
consideration are used and ^ is taken equal to 0.3^ equation (24) 
shows that the maximum damping with negative gearing corresponds to 
T 1/2 = 1.0 second. 

Figure 15 shows the curve for this value of damping. The cusp 
point for negative gearing occurs at = 3.96 radians per second and 

K = -O .035 deg^de§/sec. For positive gearing, on the other hand, an 
infinite range of cUq values which will give better than this damping 
with ^ = 0.3 is seen to be available. This bears out the previous 
statement that the design, with regard to damping ratio and natural fre- 
quency, Qf a conrpiX)mlse autopilot which will improve a variety of fli^t 
conditions is less restricted when positive gearing is used. 


Since the curve of figure 15 is typical of the curves for dampings 
somewhat greater than the airplane damping, it will be discussed in more 
detail. Thb curve has two separate branches and the discontinuity occurs 
between cu = 4.8 radians per second and cu = 5 radians per second. This 

critical frequency actually is oocr = V^o ” (see eq. (21a)). From 
what has been said previously, the portion of the curve of greatest 
interest is the second half (for positive K) . This starts at the crlti- 

P„ + 2R 

cal frequency at (Oq = «> and K =--^ This value of K ' is the 

^1 

value for which an ideal autopilot (one with ho lags) would yield the 
given damping, and will be denoted as Kq. The bottom part of the boimd- 
ary can be seen to represent the values of K and cSq which make the 
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aii^lane mode have the given damping, since the critical frequency is 
near the aiiplane frequency. On the other hand, the top part of the 
boundary represents the values of K and cdq at which the autopilot 
mode has the given dan 5 )ing, since the values of cd along this part of 
the curve are very close to the corresponding values of cUq. 

Figure 16 shows the boundary for = O.60 second. The loops 

for negative damping have disappeared. This cturve is the cusp curve 
for positive K, obtained by using equation (22) in equation (l4) . This 
fact has little practical significance, however, since there is no funda- 
menteil change in the shape of the curves at this dan 5 >ing. 

For large values of damping, further changes occur in the type of 
damping ciorves in the KniQ-plane. For small values of R, equation (21a) 
gives real solutions for oDo with any value of tn. However, for larger 
vadues of R the discriminant of equation (21a) changes sign for certain 
combinations of cn and R, so that no real solutions for oJq exist. 
Setting this discriminant equal, to zero yields 

(cu^ + R^)^ + |^^(Po + 2R)2 - 2R(Po + 2R) - 2(^ iu^ + R^) + 

Qo[^ + 2R(Po + 2R^ = 0 ( 25 ) 

p p 

When this eqiiation has a positive, real root for the quantity cd + R , 
this root determines the range of real values of cd for which the solu- 
tions for (Do in equation (21a) are complex. The value of R for 
which real roots occur in equation (25) is obtained by equating the 
discriminant of this equation to zero to give 

(1 - . (i . ^ = 0 (26) 

For the airplane being considered, with ^ = 0.3^ the values of R 
obtained from equation (26) are R = -I.58 and R = I.61. Since only 
the positive-damping boundaries eire of interest, only the negative value 
(which corresponds to T^/p = second^ need be considered. For this 

value of damping, the range of values given by equation ( 25 ) becomes a 
single value. For greater damping, equation (25) gives a finite range 
of values of cd which will not occ\rr at any real value of cDq. Thus, 
no value of cd in this range can occur for the recjuired damping when 
^ '= 0 . 3 * The constant-damping boundaries for these larger values of 
damping, which are characterized by the absence of a given range of 
CD values, are of a somewhat different type from the ccirves for lower 
damping . 
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As an example of the type of cvirve which occiars for large dan^jing 
values, the boundary for T^/2 = 0-25 second (E = -2.77) is presented 

in figure 17* Substituting 5=0-3 and R = -2.77 into equation (25) 
gives 0) = 4.3 radians per second and to = 6.3 radians per second. 

That is, no values of 0 in the range 4-3 < n> < 6.3 radians per second 
occur on this curve. Some of the significant features of this curve 
will now be described, since it is typical of the hi^er-damping curves. 

The significant portion of this curve defines the same general type 
of wedge-shaped good region as was present for lower dampings. On the 
lower part of the significant boundary, the frequency of the mode with 
the given damping is in the range 4.0 < cn < 4.3 radians per second; 
that is, the frequency varies between the critical frequency and the 
smaller frequency given by equation (25) . These frequencies clearly 
represent the aiiplane mode. At the point of ^he wedge the previously 
discussed discontiniiity appears. At this point both the airplane and 
autopilot modes have the same damping. On the upper part of the bound- 
ary the autopilot mode has, the required damping and the airplane mode 
has higher damping. 

Figure 18 is a collection of the significant portions of the damping 
curves previously discussed. This figure shows that, as the required 
damping Increases, the wedge-shaped region in the positive Kn>o -plane 
which insures this amoixnt of damping moves upward and to the ri^t. The 
figure eilso shows that the damping obtainable with a second-order yaw 
autopilot cannot always be increased merely by increasing the gearing 
ratio. The reason is that, for a given value of cDq, increasing K 

beyond a certain value makes the autopilot mode less stable. Clearly, 
more damping can be obtained by increasing cOq at the same time that 
K is increased. These results confirm the statement made in the dis- 
cussion of figure 10 that ft may be necessary to increase the natural 
frequency when the gain is increased. 

It is Important to remember that, as in the 5^J0o~Plaiie^ "two oscil- 
lations are present in the regions of most interest in the KLOQ-Plane. 

The typical frequencies given on the individual boiindaries indicate that 
the portions of the wedge-shaped boiindarles where K is high correspond 
to the autopilot mode having the given damping, whereas the parts where 
K is low correspond to the aii^plane mode having the given damping. 

Figure I8 seems to indicate that an infinite amount of damping 
might be obtainable by simultaneotisly increasing K and cuo, in contra- 
diction to the discussion concerning figures 10 and 11. However, further 
changes occur in the type of curve at larger dampings. Equation (21a) 
indicates that a change mi^t be expected when R > \/^, since the quan- 
2 2 

tity Qo - R - 0) does not go through zero for any value of co. That 
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is, there no longer is a critical frequency oicr =VQo “ Moreover, 

the good regions as previously shown do not take into account the bound- 
ary of equal roots. That is, within the good regions no oscillation has 
less damping than the indicated amount; however, a real root with less 
damping than the indicated amount may he present at points above the 
boundary of equal roots. The explanation is as follows. At the o) = 0 
point of a given hi^-damping ciorve (which point is on the bo\mdary of 
equal roots) there are two equal real roots with the indicated amoimt 
of damping. However, at points above the boundary of equal roots, one 
of these real modes will decrease in damping. The oscillatory mode 
which breaks down into two real modes at the boundary of equal roots in 
the KcDo-plane is the airplane mode, since, as can be seen from fig- 
ure 18, oscillations at autopilot frequency do occur above the boundary 
of equal roots. Thus, for points above this boundary the airplane mode 
becomes nonoscillatory. The physical reason for this phenomenon is that 
large gearing ratios ca^lse the airplane mode to be overdamped. The 
boundary of equal roots in the ^cnQ-plane, on the other hand, corresponds 
to large ^ values. Therefore, as can be seen from the curves, the 
mode which is overdamped at points above this boiondaiy is the autopilot 
mode . 


Strictly speaking, the boxandary of equal roots should be considered 
as the upper limit of the good regions shown in figure I8. Actually, 
this restriction is necessary only for the regions of very hi^ damping. 


For example, when the point K = 0 . 


60 degyiieg/sec 


and U)q = 21.5 radians 


per second is taken on the boundary for T]^/2 = O.60 second in fig- 
ure 18 and the characteristic roots are found, the autopilot mode has 
the given damping ~ 0 > 6 o second at o> = 21 radians per second, 

whereas the airplane oscillation breaks down into two well-damped non- 
oscillatory modes with '^^ 2-/2 ~ 0*22 second and Tj ^/2 ~ 0.09 second. 

Thus, for this moderate damping, points above the boxandary of equal 
roots still give the required damping. 


For extremely large values of dan^Jing, however, the boundary of 
eqioal roots becomes important. In fact, the maximum damping of the 
complete system at a fixed value of ^ is the largest negative value 
of R occurring on the boundary of equal roots. The points on the 
botmdary of eqvial roots are obtained by tising o> = 0 in equations (21a) 
and ( 21 b) . The largest vadue of R occurring on the boundary is 
obtained by setting the discriminant of equation (21a) eqxaal to zero 
with oi = 0 . This procedure gives the quartlc in R (from eq. ( 25 ) 
with £0 = 0) : 


(4^^ - 3 )r^ + 2Po(2^^ - i)r 3 + + 2Qo)r^ + 2P0Q0R + Qo^ = 0 (2?) 



MCA TU 2857 


25 


The largest negative real root of this equation is the value of the 
maximum damping. For the example considered, this value is E = - 5 .11 
and the corresponding T^/g ~ 0 .l 4 second. For this value of damping, 
the good region becomes a single point on the boundary of equal roots. 

The airplane mode is a double real root equal to - 5 . 11 ^ and the auto- 
pilot mode is an oscillation with the same amount of damping. The 
autopilot characteristics for this maximum- damping point were used in 
calculating the motion subsequent to a 5 ° sideslip disturbance (fig. 19 ). 
For practical puaprases, the aii^lane motion can be considered as non- 
oscillatory, since the autopilot oscillator is of too hi^ a frequency 
to have an appreciable effect on the aiip)lane. The well-damped nature 
of the Dutch roll motion (which has now become nonoscillatory) is evi- 
dent in the 3 -motion. The slow return of the roll motion is due to the 
spiral mode. Although the li^tly damped spiral mode is generally not ’ 
considered troublesome, it is necessary to keep in mind that the dis- 
cussion in this paper deals only with the improvement of the damping of , 
the Dutch roll mode, and that the two aperiodic modes in the airplane’s 
lateral motion have been ignored. 

Effect of varying ^ on curves in the KiDq - plane.- In order to 
obtain an idea of the effect of changing ^ on the curves in the 
Kuo-plane, a comparison of the boundary of equal roots, the alip>lane 
damping boundary, and the boundary for T-j_/2 = O.60 second is presented 

in figure 20 for ^ = 0 . 3 ^ 0 . 6 , and 0 . 9 * Althou^ the region of improved 
damping with negative gearing increases in size, the narrowness of the 
frequency range and the other difficulties previously mentioned still 
make the use of negative gearing impractical. The increased slope of 
the upper pajrt of the boundaries siii5)ly implies, as would be expected, 
that for larger values of leirger values of K are required at any 
value of a>o to make the autopilot mode become unstable. The variation 
in the position of the boimdary of equal roots simply implies that the 
more highly dan5>ed mode associated with each point on a given boundary 
tends to become critically damped at lower gearings as ^ increases. 

None of these variations appear to be very important practically. 
Indeed, the most important fact about the effect of varying ^ is that 
the bottom part of the boundary approaches the same value of K at 
laorge values of o^o for a 1 1 valuesi of This behavior is due to the 

fact that the asyn5)totic value of K is Kq, the value required to 
obtain the given amount of damping with an ideal autopilot. Now, the 
bottom part of the boundary can be seen to be the in^jortant part, since 
it gives the lowest value of K for which the system has the reqiiired 
damping; also, the mode which attains the required damping at the points 
on this part of the boundary is the airplane mode. The flatness and 
invariance with ^ of the bottom part of the boundary imply that the 
mlnimimi value of K for which a given damping may be attained is rela- 
tively invariant for changes in ^ and cDq at fairly large values 
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of ajQ. This fact will he iisecL in the next section to obtain a simple 
rule-of-thumb method for designing an efficient yaw-rate damper which 
will insxire a required (damp ing for several extreme fli^t conditions of 
an airplane using low gearing. 


Design of Compromise Yaw-Eate Damper 

For the sake of simplicity, the autopilot to be designed will be 
assumed to provide only a yawing moment proportional to yawing velocity, 
as in the previous analysis. It may therefore be called a 0^^ auto- 
pilot. Values of K, ancb, (Dq which will efficiently improve the 

damping of an aiip»lane in various fli^t conditions are desired. The 
criterion of the efficiency of the autopilot will be that the gearing 
ratio required for the autopilot must be small. This means that the 
autopilot power required will be small. Moreover, the use of small 
gearing makes the yaw damper easier to override in steady turns. 

Table II gives the parameters used for the three fli^t conditions 
of the airplane chosen for the present exang)le. Case I is a hi^-lift- 
coefficient, low-wing- loading landing condition at sea level. Case II 
is a low-lift-coefficient, medium-wing-loading cruising condition at 
30,000 feet. To con 5 )lete the plctxire, case III is a hi^-lift- 
coefficient, hi^-wing- loading cruising condition at 30,000 feet. These 
cases will serve as examples to illustrate the metho<3. 

Table II shows that the Dutch roll oscillation in case II is very 
poorly damped, since it requires almost 7 seconds to damp to one-half 
amplitude. Although the other two cases are not so poorly damped, they 
are still unsatisfactory. Calcxilated motions for the three cases in 
response to a 5° disturbance in sideslip are shown in figure 21. No 
attempt will be made to set up any con^plicated criteria for adequate 
damping. Instead, the criterion chosen for purposes of illustration 
will be that the Dutch roll oscillation should damp to one-half ampli- 
tude in 1 second or less at any flight condition. Actually, the auto- 
pilot may be designed to insure a different amount of damping for each 
fli^t condition, in case one of the flight conditions is required to 
be more stable than another. 

The fimdamental problem is to find the set of points (K,oio, which 
satisfy the given damping criterion for all three fli^t conditions. 
Moreover, the mlnimimi possible gearing is desired. The value of Kq 
for Tj^y'2 “ ^ second is calciilated for each flight condition. From 
the eqixivalent-oscillator analysis, the value of Kq is shown to be 


Ko = 


1.386 
^1 '^1/2 


- Po 
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Substitution from table II gives the values of Kq as 0.1213, 0.0698, 
and 0.137^ deg^deg/sec for cases I, II, and HI, respectively. 

If the largest of these vadues of Kq (corresponding to case III) 
is chosen for the geeiring of the autopilot, the curve for T1/2 = 1 second 
in the ^oio-plane will be the critical curve for case III, but for cases I 
and II the curve for Tj/2 ~ ^ second will be a curve between the air- 
plane damping and the critical damping. Each of these three c\arves 
defines an infinite region in the upper right portion of the ^a>Q-plane 
as the region insuring better damping. Therefore, in the -plane an 
infinite region of points which are common to these three regions must 
exist and can be used with this value of K to obtain better damping 
than T]_y2 “ ^ second for all three flight conditions. A plot of these 

regions is shown in figure 22(a), where K = 0.l4 deg/deg/sec has been 
used. The plot shows that any set of values in the -plane which 
insures ~ ^ second for case II, the hlgd^-fnequency case, will 

also insure this amount of damping for the other two cases. However, 
any point in this region is at a value of cdq considerably above the 
optimimi for cases I and III. As was previously pointed out, for such 
values of (Oq the TnlninniTn value of K for which a given damping may 
be obtained does not differ much from Kq. The Implication is that the 
minimum value of K for which ^ second may be insured for 

case III, while at the same time T1/2 remains less than 1 second for 

case II, is not much less than 0.1374 deg^deg/sec. To confirm this 

hypothesis, the necessary points on the curves in the 

Tj^^ 2 “ 1 second were obtained for K = 0.12 deg^deg/sec and 

K = 0.10 de^deg/sec, and the significant parts of these curves are 
shown in figures 22(b) and 22(c). 

For values of K < 0.1374 deg^deg/sec the curve for T-j^^2 = ^ second 
for case III must be a loop, since the damping is higher than the critical 
damping. Figure 22(b) shows that when K = 0.12 de^deg/sec this loop 
is still large enough to intersect the region where ^ second for 

case II. The hatched region on figure 22(b) is the region which will 
ins\ire ^ second for all three flight conditions at K = 0.12 

deg/deg/sec. Figure 22(c) shows that the loop for case III has become 
too^ small to Intersect the region where ^ second for case II 

when K = 0.10 deg/deg/sec. Therefore, a value of ^ second 
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cannot "be obtained for cases H and III simultaneously when K = 0.10 
deg/deg/sec. The c\irve for case I is not shown in this figure since, 
as seen in figures 22(a) and 22(h), it is not necessary. 


The discussion of figure 22 shows that, in order to improve all 
the f 1 1 ght. conditions in the present example with a single second-order 
autopilot, a gearing ratio has to he used which is almost as large as 
would he needed for an ideal autopilot that would stabilize all the 
f 1 1 ght. conditions. This difficulty arises because the low-frequency 
conditions require the hi gb est values of Kq, as can be seen by 
examining the approximate expression for the frequency of the Dutch 
roll oscillation at any fli^t condition and the associated expression 
for Kq: 



^^o = 


- ^)(P - Po) 

% 

Cng 


fa 


(P - Pq) 


( 28 a) 


( 28 b) 


Becaiise of the possibility of variations in 0^^, and P - Pq 

in the various flight conditions, it is not possible to state that in 
general the low-frequency fll^t conditions will require the hl^er 
values of Kq, as is true in the present example. When the highest 
required value of Kq occurs at the hi^-frequency conditions, the 
Tnlnlmiim value of K for the compromise autopilot may be considerably 
below this hipest required value, because the region of overlap similair 
to that shown in figure 22(b) will include the optimum point for the 
high-frequency condition. 

In any case, the characteristics of the compromise autopilot can 
be obtained by plotting the regions insuring the required damping for 
the various conditions at various veilues of K below the maximum value 
of Kq, as is done in figure 22. The value of K which leads to a 
small overlap region is then the mlnimimi compromise gearing, and the 
values in the overlap region of the ^o>Q-plane define the possible values 
for ^ and cDq of the compromise autopilot. 



MCA TN 2857 


29 


A convenient method is available for detennining the range of 
values of K to he used in obtaining the overlap region. The veilues of 
Ko and the optimum value of K which will yield the required damping 
for each of the flight conditions aire calculated. Then the con^promise 
value of K must lie between the hipest optimum value of K and the 
highest value of Kq. In making the plots it may soon become clear that 
some of the flight conditions do not need to be considered, as happened 
with case I .in the present exanqjle. 


In the present example, it has been shown that a gain much smaller 
than the largest value of Kq necessary cannot be used, unless a vari- 
able gain is available. If a constant gain must be used at all flight 
conditions of such an airplane, a simple rule of thumb for designing 
the autopilot would be to choose the largest value of Kq (corre- 
sponding to the low-frequency conditions) and values of ^ and coq 
near the optimum point of .the hi^-frequency condition. In this way 
some advantage is derived from the second-order characteristics of the 
autopilot in that the damping obtained from the autopilot is much larger 
for the hi^-frequency conditions than that which would result from an 
ideal autopilot with this gearing. 


In applying this sinple method, 0 


.l4 deg/deg/ 


sec was chosen for 


the value of Kj and the optimimi point was fotmd for case II . The opti- 
mtom point was ^ = 0.523 and = 9*^9 radians per second. Figure 23 
shows the calculated motions for the three flight conditions with this 
autopilot, subsequent to a 5° sideslip disturbance. A comparison of 
figure 23 with figure 21 reveals that the stability of all three flight 
conditions is improved. In particular, case II is greatly improved, 
while cases I and III both have T-i /p sli^tly less than 1 second. 


Since the results of figure 21(b) show that a gearing of 
0.12 de^deg/sec could have been used, the optimum point for case II 
was calculated for this gearing and found to be ^ = 0.k8^ and 


03r 


= 8.81 radians per second. Figure 2k shows the motions in the three 


flight conditions with these autopilot characteristics. A comparison 
of this figure with figure 23 reveals that the danping obtained with 
this autopilot is only slightly less than the damping obtained when 

K = O.l^J- deg/deg/ sec. 


- w/ - 


If the low-frequency conations require little improvement, the 
problem approaches that of improving only a single fli^t condition, 
so that the optimum-point characteristics may be used to decre.ase the 
gearing necessary. For example, suppose that Tj/p = 1.5 seconds had 


The 


been considered satisfactory for the low-frequency conditions, 
value of Kq necessary to obtain this danplng is 0.075 deg^eg/sec. 
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and the corresponding optimum point for case II is ^ = 0.389 and 
(Dq = 7-^ radians per second. The improvement of the stability of the 
various fli^t conditions may he seen from the motions presented in fig- 
ure 25 . In these motions the value of T]^/2 cases I and III is 

between 1 and 1.5 seconds, hut the value of T^/g case II is well 

under 1 second. 

It is cleax from the previous discussion that for the present 
example an autopilot with fixed characteristics would require a leirge 
gain (corresponding to the largest value of Kq) in order to improve 
both 1 anding and cruising fli^t conditions. A variable gain would 
therefore be desirable, so that this excessive gearing would not have 
to be used at the hi gh er speeds . Althou^ it would be impractical to 
expect that all three autopilot parameters should be variable in fli^t, 
as woiild be desirable to obtain the optimum autopilot for each fll^t 
condition, it would be relatively single to make the gain variable. The 
value of 5 chosen would be that of the optimum point for the high- 
frequency condition (as calculated for the low-gearing value), and the 
value of 03o would be preferably sli^tly above the optimum value of 

cDq for the high-frequency condition. (Because of the rapid change in 
damping at values of oiQ less than the optimum, in practice a safety 
factor should generally be added to this value.) In such cases, prob- 
ably only two gain positions would be necessary - a hl^ gain for the 
low-velocity conditions and a low gain for the hig^- velocity conditions. 
Thus, low gain coTild be used in the hi^- velocity flight conditions, so 
that the adverse effects of the autopilot in steady turns at high 
velocity would be small. 


Validity of As sinning Pure Cjj Autopilot 


A few final remarks will be made concerning the assumption that the 
autopilot is sensitive to yawing velocities only, and provides yawing 
moments only. For practical autopilot installations , the sensing device 
(usually a gyroscope) is fixed in the airplane. The device is therefore 
sensitive to yawing velocities about some axis fixed in the aiiplane . 

The equations of motion, however, are set up with respect to the sta- 
bility axes. For various flight conditions, the angle of inclination 
between the gyro axes and the stability axes varies because of the 
varying angle of attack. This angle of inclination makes the gyroscope 
sensitive also to rolling velocity about the stability axes. In addi- 
tion, the displacement of the autopilot-actuated control surface from 
the longitudinal stability axis gives rise to rolling moments. As shown 
in reference 10, these two effects cause Increments in 

C 2 due to the autopilot, in addition to the expected 


‘P 
The 


and 
Increment . 




effect is the most important of these in affecting the Dutch 


roll stability of the airplane, and this effect will now be discussed. 
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When the stahllity of a single flight condition is to he improved, 
the Cjip effect may he removed hy orienting the gyro axis along the 

flight path, and no modification of the results is necessary. If an 
aingle of gyro inclination does exist, the general method of setting up 
the constant-damping boundaries is, of course, still valid. In order 
to take account of the effect, Dijr is simply replaced in equa- 

tion (6) hy the proper linear combination of DiJr and T^, as shown in 
reference 10. Because of the presence of the D(^ term all three degrees 
of freedom must be used, and the equivalent-oscillator simplification can 
no longer be used. Thus, the method is more complicated when the 

Sr 

effect of the autopilot is not small enough to be neglected. When the 
Cnp effect becomes excessively large, the problem becomes even more 

complicated because of the fact that the damping of the aperiodic modes 
becomes important, and curves of constant damping for these modes (corre- 
sponding to real roots) have to be plotted. In fact, the aperiodic modes 
may combine to form another oscillation. 


In considering the problem of simultaneously stabilizing various 
flight conditions which have varioTis angles of attack, the Cn^ effect 

must always be considered unless the gyro axis can somehow be rotated 
so that it is always parallel to the relative-wind axis. The effect 

can sometimes be ignored, if the angle -of- attack range of the fli^t con- 
ditions is small, by choosing an orientiation of the gyroscope in such a 
direction that all the angles of inclination are small. Because of the 
extreme conplexity of the problem when the C^p effects must be con- 
sidered, this paper is confined to the consideration of Cn^. effects 
only. 


COrfCLUDIWG KE 3 MRKB 


The damping of an oscillatory system that makes use of a second- 
order rate damper with a given geeiring ratio K can be improved by 
adjusting the shape of the autopilot frequency response (that is, the 
damping ratio and the natural frequency of the autopilot) . For the 
piirpose of determining the effect of a second-order yaw-rate damper on 
the damping of the Dutch roll motion of an airplane in a given fli^t 
condition, the aiiplane may be represented as an equivalent oscillator. 

By using this equivalent-oscillator concept, the optimum shape of the 
autopilot frequency response corresponding to a given gearing ratio or 
required damping may be obtained from a simple set of equations for any 
fli^t condition. The geeiring ratio necessary to obtain a given amount 
of damping when the damping ratio and natural frequency of the autopilot 
are near their optimum values is considerably less than the gearing ratio 
necessary to obtain the same amount of damping with an idealized (no-lag) 
autopilot . 
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The prohlem of designing a second-order yaw damper for an airplane 
which requires improvement in damping for several fli^t conditions is 
more complicated, since for each fli^t condition the airpleine is repre- 
sented by 'a different equivalent oscillator. However, a simple method 
of compromise is derived for fli^t conditions in which the effect of 
the yaw-damper sensitivity to rate of roll is small enou^ to be 
neglected. 

The effects of a second-order yaw damper on the stability of any 
given fli^t condition of an airplane can be obtained by examining the 
constant-damping curves in the plane of the damping ratio and natural 
frequency of the autopilot and in the plane of the gearing ratio and 
natural frequency of the autopilot. Theoretically, any given fli^t 
condition may be stabilized by using either positive or negative gearing. 
When negative gearing is used, the autopilot natural frequency must be 
less than the aiiplane frequency. The use of negative gearing is shown 
to be impractical, however. For positive gearing, the autopilot natural 
frequency must be greater than the airplane frequency. 

For fixed positive gearing, there is an infinite number of combi- 
nations of autopilot natural frequency and damping ratio for which the 
second-oitler autopilot gives better damping than an ideal autopilot of 
the same gearing. For. fixed positive danping ratio of the autopilot, 
there is a range of values of positive gearing ratio which will provide 
a given damping to the system at any autopilot natural frequency greater 
than the airplane frequency. Increasing the gearing ratio of the auto- 
pilot to excessive values will always cause the autopilot mode of oscil- 
lation to become vinstable for a given damping ratio and natural frequency 
of the autopilot. If larger gearing ratios are to be used in order to 
obtain hi^er stability for the airplane mode, the damping ratio or the 
natural frequency of the autopilot, or both, must be Increased. Expres- 
sions are derived for the maxiimnn damping' under various conditions . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 2, 1952. 
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APPENDIX A 

EXPRESSIONS FOR CONSTANT -DAMPING CURVES AND DISCUSSION 
OF EQUIVALENT-OSCILLATOR CONCEPT 


Since the characteristic eqioation of the aii^lane-damper system is 
of sixth degree, as given in equation (7)^ i"t may he written 

F(d) = D^ + A^D^ + A^D^ + A 3 D 3 + A 2 D^ + A^D + Aq = 0 (Al) 

Here A^, A]^, . . . Aq are functions of K, and for a flight 

condition in which the airplane parameters in equations (3) to (5) are 
known. If any one of the three autopilot parameters is fixed, stahillty 
boundaries may he obtained in the plane defined by the other two param- 
eters by letting D = R + icu, fixing the value of R for each curve, 
and varying oi. Since equation (Al) is a complex equation, it may be 
written as two real equations and solved for the two autopilot variable 
parameters at each value of D = R + Ixii. 

According to reference 6, if 


fn(H) = ^ , (n = 1,2,3,...) (A2) 

dR^ 

then the two real equations obtained from equation (Al) by setting 
D equal to R + l £0 are 

F(R) - 03^f2(R) + u>^f4(R) - CD^ = 0 (A3) 

f^(R) - co2f2(R) + cD^f5(R) = 0 (A4) 

Since the coefficients of • F(R) and fn(H) are functions of the 
two variable autopilot parameters only, the values of these two param- 
eters which will yield an oscillatory mode of motion with a given 
damping may be obtained by choosing a value for R and solving equa- 
tions (A3) and (a 4) at any value of u). This procedure yields a point 
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on the given E curve in the parameter plane; for instance, in the 
^cjOo -plane if ? and coq "be chosen as the variable parameters. In this 
case, calculations are simplified hy the use in equations (A 3 ) and (a 4 ) 
of ,A and B as defined in equation (l), since these equations become 
linear in A and B. 

Because of the increasing difficulty of obtaining accurate values 
of stability derivatives for present-day airplanes, the characteristics 
of the airplane may be given in terms of transfer functions which are 
obtained from fll^t data. In order to obtain the form of the char- 
acteristic equation in terms of transfer fmctions, equation ( 7 ) is 

divided throu^ by (d^ + 25o5qD + 05 q^)Fq(d) . The resulting equation is 




D^ + 25 cDoD+cDo^ 


The first factor in the second term of this equation is the autopilot 


transfer function 


li’tjP 


(d) (See eq. 8 .) If equations (3) to ( 5 ) are 


solved operationally for zero initial conditions, the expression in 
brackets is foimd to be the airplane transfer function: 


Dt 


UA 


(D) 


— - C>*i« — — — ^ — ■■■■■■■ I ■ - . ■ - ^ ^ ~ ~ — I— ■ 


2 


Fo(D) 


(A5) 


Therefore, the characteristic ^equation may be written 


1 - 




= 0 , 


(a6) 


in terms of the autopilot and airplane transfer functions. The airplane 
transfer function for the lateral motion, as shown in equation (A 5 ) , 
should have -the form of a cubic in D over a quartic in D. Therefore, 
when equation (a 6 ) is cleared of fractions a sixth-degree equation of 
the form of equation (Al) is again obtained as the characteristic 
eqtiation . 



MCA TN 2857 


35 


The setup and solution of eqiiations (A3) and (a 4) woiild be greatly 
simplified if the characteristics of the airplane lateral motion coiald 
be expressed as a second-order system rather than a fourth-order system. 
Since the yaw damper is used essentially for damping the Dutch roll 
oscillation, it was thou^t desirable to see whether the airplane could 
be replaced by an equivalent oscillator whose damping and period were 
those of the Dutch roll oscillation of the airplane. That is, if the 

Dutch roll mode, as calculated when three- degrees of freedom are con- 

o 

sidered, is obtained from the quadratic equation D + P^D + = 0, 

the question is whether the aiiplane characteristics may be represented 
by the oscillator described by 




2 %Z U^2 


K, 


2 i<P 


o) = 


Equation (A5) may be rewritten in the form 


Di^ 


L-^-lA 


(D) = 




1 

0 

60 

(D + a^) 

(d^ + a^D ^ 

^ 3 ) 


C 2 ^XZ^'^ 
^ ■ Kx^ y 

i(d2 + 

D^D^ 

+ + a^] 



-Cj^D 

^D + a^^ ^D^ + + a^^ 

d2 + p^d + 

D^D^ + aijD + 


(A7) 


where the constant is defined by 


^ 

V % 


The expression in the first bracket on the right-hand side of equa- 


tion (A7) is the same as woiild be obtained for 



if the assump- 


tion were made that the ai2rplane could be represented by the equivalent 
oscillator described pre-vLously. Thus, if the expression in the second 
bracket has a frequency response which does not significantly differ 


o 
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from unity for the range of ifrequencies which is- important in the air- 
plane frequency response, then the expression in the first bracket cem 
he used a,s a valid appix)ximation to the aiiplane frequency response. 
Therefore, the expression in the second bracket can be considered as a 
correction fmction. 


The factor + a]^D + a^ yields the two real characteristic roots 

of the airplane lateral motion corresponding to the spiral mode and the 
damping- in-roll mode. The us\ial approximation to the damping- in-roll. 


root is 


VC2p/4bp^/, 


and the spiral root is of the order of 


VCYp/2bp^. 


However, the numerator 


of the second bracket can be shown to be 


d3 


v/^P 


+ 




y3 


If Cjp = 0, this cubic has the roots 

Cy 


D = 0 


2H^ b 


D = 


V 


which are merely approximations to the real roots which occur in the 
cubic in the denominator of equation (A7) • li can therefore be expected 
that the second bracket will usually be close to unity and that no large 
errors will arise if the airplane lateral motion is represented by a 
single-degree-of -freedom oscillator in yaw with nbndimens tonal inertia 


same period amd damping as the Dutch roll oscil- 
lation. This is especially time for small values of 


For the airplane used as an example in this paper, the correction 
. (B + an)(D^ + apD + ao) 

function ^ was evaluated throu^out the fre- 

D^D^ + aj^D + a^^ 

quency range and its effect on the total airplane transfer function was 
found to be very small except close to zero frequency. Finally, actual 
comparisons of the stability boundaries for the three-degree-of-freedom 
aneilysis and the equiveilent-osclllator approximation showed that the 
approximation was valid. 
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When eqiiation (A7) ±s used with this approximation, the 'character- 
istic equation as given in equation (a 6) becomes 


1 + 


04 

^0 

1 

+ S^oioD -1- o>Q^ 

•f 

+ 

04 

-a_j 


= 0 


or 




(A8) 


The characteristic equation of the -system is therefore a quart ic, 

F(D) = + (Pq + 2^03o)d3 + (Qo + -f 2 ^o>oPo)d 2 -f 

s -t- (Pq + a)d^ + (Qq + B + APo)d^ + 

(PqB -1- QqA -t- CjKB^D + QqB 


= 0 


(A9) 


The equations for the stability boundaries which replace equa- 
tion (A3) and equation (aA) for a fo\jrth-order characteristic equation 
are 


F(K) - + (JD^ i= 0 

fl(R) - o^f3(E) = 0 

The functions fn(R) are obtained from equation (A2) 


(AlO) 


(All) 
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As an example, tlie actual parametric equations used for the equiva- 
lent oscillator will he derived. Equation (A9) is used in equation (A2) 
to give 

fl(n) _ = 4 r 3 + 3R^(Po + A) + 2R(Qo + B + APq) + (Pq + CiK)B + Q^A 
fgC-R) = 6 r^ + 3R(Pq + A) + Qq + PqA + B 
f^(R) = 4 r + Pq + A 


These expressions are used in equations (AlO) and (All) . To obtain the 
constant-K curves (that is, curves in the AB-plane or ^oJq - plane), the 
resulting equations are written in A and B, as follows: 


[(3R + Po)o:^ - r(r^ + PoR + Qo^A + - [r2 + (Pq + CiK)R + 

03^ - (6r2 + 3PqR + Qq)o3^ + r2(r2 + p^ + Qq) 


B = 


(A12) 


- (3R^ + 2PqR + Qq^A - (2R + Pq + CiK)b = 
-(4R + Po)a3^ + r(^2 + 3PqR + 2Qo) 


(AI 3 ) 


For a given value of K, any constant-damping curve is obtained by 
fixing R at the appropriate value and taking a sufficient nuBiber of 
positive values for a> to obtain the points necessary to determine the 
cxirve. The resulting two linear equations in A and B are solved to 
obtain each point. If cxirves in -the ^cOq - plane are desired instead, the 
values of ^ and ojo can be obtained for each set of values of A 
and B for which B is positive. 

The constant-^ curves in the Ko>Q-plane are obtained by solving the, 
parametric equations for K and cdq instead of for A and B. For 
exanq)le, solving equation (AI 3 ) for K gives 
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K 


" + (Qo + 2PoR + 3 R^)a + 

R(2Qo + 3PoR + 4r2^J^ 

= — ^ -s (2 ?cdo + 1*0 + 4r)cd 2 - F(Po + 2R)05 o^ + 2^0)o(qo + PqR + 3R^) + 

L L 


r(2Qo + 3PqR + 


(A14) 


and using this eqxiation in eq-uation (A12) gives 


(Qo -R^- o)^)b - (Po + 2 R)(a )2 + R^)A+ai^- (Qq + 2 PqR + 2 R^)(a^ + R^) = 

(Qq - R^ - - 2^0)3 (Pq + 2r) (u>^ + R^) + co^ - (q^ + 2PqR + 2R^)(a>^ + R^) = 0 

(AI5) 

This eqiiation is a quadratic in cdq for fixed values of ^ , R, 
and 03. The values of 03o from equation (AI 5 ) are then used in equa- 
tion (A14) to obtain the points on the curves in the K03o“Pl^®* Since 
equations (A12) and (AI 3 ) are linear, there is only one point in the 
AB-plane or the ^03Q-plane for, each oscillatory mode when K is fixed. 

However, for a fixed value of ^ there may he two points in the 
KtDo-plane for a given oscillatory mode when the quadratic equation (AI 5 ) 
has two reail, positive roots for' (Dq. • 
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APPENDIX B 

DEEIVATION OF EXPRESSIONS FOR MAXIMUM-DAMPING, 
CHARACTERISTICS WITH FIXED GEARING RATIO 
BY USE OF EQUIVAIENT-OSCILLATOR CONCEPT 


The characteristic equation of the airplane — yaw-damper system, 
for the assun 5 >tion that the airplane may he replaced, hy an equivalent 
oscillator, is 

D^ + (Pq + a)d3 + (Qq + b + APo)d^ + (PqB + AQo + CiKb)d.+ QqB = 0 

(Bla) 

The conditions on A and B which result in the maximum damping for 
a given gearing ratio can he calculated as follows. For the case of 
maximimi damping, equation (Bla) will have two pairs of equal complex 
rootsj hence this equation becomes: 

(d^ + PD + q)^ = D^ + 2PD^ + (p^ + 2 q)d^ + 2PQD + (Blh) 

o 

where D + PD + Q is the quadratic corresponding to the double oscil- 
latory root at the cusp point. The following four equations are obtained 
by equating like coefficients in equation (Bib) : 


2P = Pq + A 

(B2) 

P^ + 2Q = Qq+B + APq . 

(B3) 

2PQ = PqB + AQo + CiKB 

(b4) 

= QqB 

(B5) 


Since K, C^^ Pq^ an<i Qq are known, these four independent equations 

may be solved for A, B, P, and Q. Thus, the autopilot parameters 
which will give the hipest damping to the Dutch roll oscillation, and 
also the period and dainping of this oscillation, may be determined. 
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From eqmtions (B2) and (B5) 

I 


A = 2P - Pq 

(b6) 


(B7) 

Qo 


Using equations (b 6) and (B?) in equation (B4) yields 


(Pq Cj^K)Q^ - 

2«o(9 - 9o) , 


(B8) 


When equations (b 6) and (B?) are used in equation (B3 ) , the result may 
he written 


9c(P - ^of ' (9 - 9o)^ 


(B9) 


and hy stibstituting equation (b 8) into equation (B9) a quartic equation 
in Q is obtained, which may be written 




Fo(9 - 9o)^ + 


= H,(9 - 9o) 


Finally, this equation yields the two quadratics in Q: 


'1 + 


C]_K 




- 2QqQ + = 0 


(BIO) 


Evaluation of the discriminant of equation (BIO) reveals that equa- 
tion (bio) has only two real raots, one less than Qq and the other 

greater than Q^. As can be seen from figure 6, the smaller value of Q 

(which yields a smaller valiie of ojo) corresponds to the breakdown of the 
unstable loops fpr positive K. The larger real root of equation (BIO) 
is therefore the one which yields the maximum- damping point for positive 
gearing ratio. Negative gearing ratios will be discussed subsequently. 

The corresponding vadue of P is obtained by substituting this 
larger root into equation (B8) , and the maximum- damping point in the 
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AB-plane is easily obtained from equations (b 6) and (B7) • Thus, the 
value of the TnnYinmm damping of the Dutch roll oscillation obtainable 
with a given gearing ratio, and also the vali;es of autopilot frequency 
and damping which will yield this maximum damping, may be obtained by 
simply solving a quadratic equation (eq. (BIO)). 

The procedure is even simpler when the required damping is given 
and the values for the autopilot parameters which would make this 
damping the maY-lmnTti damping are to be determined. The value of P is 
determined by the required damping, since 


P = -2R 


1.386 

^1/2 


(Bll) 


and the corresponding optimum value of A may be obtained from equa- 
tion (b 6) . Equation (B9) now gives 

Q = Qo i \/^ (I’ - Po) (B12) 

Both vajLues of Q obtained from eqtiation (B12) correspond to optimum 
points, since they both coirespond to the given positive value of 
damping. These vadues may now be used in equation (BT) to obtain the 
corresponding valioes of B. Finally, the two values of gearing ratio 
may be obtained from equation (b 8) , which gives 


K = 


2PQ - P^ - AQq 


C]B 


(BI3) 


The smaJLler values of Q (obtained from using the minus sign in 
eq. (B12)) res\ilt in negative values of gearing ratio and values of cQq 
less than the aiiplane frequency. The possibility of using negative 
gearing to improve the damping is discussed further in the body of the 
paper. 
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TABLE I 

PARAMETERS OF AIRPLANE USED AS EXAMPLE TO ILLUSTRATE 
USE OF CONSTANT-DAMPING CURVES 

Altitude, ft 30,000 

Wing loading, Ib/sq ft 65 

V, ft/sec 797 

L, ft 28 

Cl 0.23 

v-1) 80.7 

Kx^ 0.0097 

0.0513 

Kxz -0.00145 

C^p, per radian -0.40 

Cjj., per radian O.80 

Cnpj per radian -0.02 

Cj^^, per radian -O.hO 

Cy , per radian -1.0 

p 

Cnp> per radian 0.25 

Cjp, per radian -0.126 

Cng, per radian . . , -O.163 

Period, sec I.3 

T1/2, sec 2.6 

Po^ sec"l 0.537 

Qq, sec"^ '. . 23.84 

C]^, sec"^ 15.98 
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TABLE II 

PARAMETEES USED FOR THREE TYPICAL FLIGHT CONDITIONS OF THE 
AIRPLANE IN DESKm OF COMPROMISE YAW-RATE DAMPER 



' 

Case I 

Case H 

Case HI 

Altitude, ft 


Sea level 

30,000 

30,000 

Mach number 


0.22 

0.75 

0.5 

Wing loading, Ib/sq ft . . 


• 

65 

85 

^^b 


25.2 

81.25 

106.3 



• • • 

0.0081 

.0.0069 . 

0.0051 

• 


0.0433 

o.o 4 i 9 

o.o 4 o 9 

%z 

« * • 

0.0027 

0.0025 

0.0010 

Cl 


0.765 

0.262 

0.771 



• « • 

0.205 

0.205 

0.212 

% 

• * • 

-0.099 

-0.107 

-0.095 



• « • 

-0.930 

-0.878 

- 0.884 

C2 

ip 

• • « 

- 0.425 

-0.474 

-0.435 

C2 i 

• • • 

0.288 

0.200 

0.300 

% • 

• • * 

0.003 

0.010 

0.003 



• • « 

-0.165 

- 0.150 

-0.165 

^“6 * 


-0.163 

-0.163 

-0.163 

Period, sec 


2.3 

1.4 

' 2.1 

T^/2, sec 

• • k 

2.0 

6.9 

2.4 



• • • 

0.704 

0.200 

0.573 

Qq, sec"^ 

• • • 

7.79 

■ 21.4 - 

8.78 

C]^, sec“^ 

• • • 

5:63 

17.0 

5.92 

V/b, sec“^ 


8.68 

26.64 

17.77 




a 


Figure 1 .- Values of autopilot parameters A and B for which a neutrally- 
damped oscillation will occur in the lateral motion. Three degrees of 

lateral freedom; K = O.086 deg/deg/sec. Values of m in radians per 
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Figure 2 .- Values of autopilot parameters A and B for vblcb an oscil- 
lation will occur in the lateral motion with damping equal to that of 
the natural Dutch roll oscillation. Three degrees of lateral freedom; 

K = 0.086 deg/deg/sec. Values of ce in radians per second are given 

at representative points. The dashed line indicates the value of B 
copx*s Spending "to sirplsiis fi^sQucDcy# 
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I Figure 3-- Values of autopilot parameters A' and B for wMch. an oscil— 

j _ lation vill occur in the lateral motion vlth T^/2 * ^*3 secondB, a 

; typical damping value between the airplane damping and the critical 

damping. Three degrees of lateral freedom} K i* O.O86 deg/deg/sec. 

I Values of 01 in radians per second are given at representative points. 

I 

I 
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rXKLU’C *t*- VOJ.W.^D ui.ukww^---uw w 

latlon will occur in the lateral motion with T]_/2 = 0-75 second, the 
— fTh-ppp flecrreea of lateral freedom; K = O.086 deg/deg/sec. 

Values of 03 in radians per second are given at representative points. ^ 
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Figure 5.- Values of autopilot parameters A and B for which an oscil- 
lation will occur in the lateral motion with = O.6O second, a 

typical damping curve for greater damping than the critical. Three 
degrees of lateral freedomj K = O.O86 deg/deg/sec. Values of co in 
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u-t ^ Cl's* j. 0^11 (.fa. u V p's..* xu od , 


WACiA TN 2857 



B.O.E.R. 



MCA TN 2857 



COo, fod tans /see. 

Figurs 7'“ Comparison of negions in tli6 ^oiQ-piana t-liat- insure various 

amoimts of damping as calculated "by the three -degree-of -freedom analysis 
and the equivalent -oscillator analysts. K = O.O86 deg/deg/sec. 
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K o 0.086 deg/deg/secj 
= 10.66 radians per 
second; ^ = 0.11. 


(b) K = 0.086 deg/deg/sec; 
ujo = 10.66 radians per 
second; ? = 0.19^5- 


(c) K = 0.086 de^deg/sec; 
= 10.66 radians per 
second; ^ = O.503. 


Figure 8. - Comparison of motions vith three different autopilots whose characteristics lie 
on a curve for ^1/2-“ ^.6o second (corresponding to points 1^ 2_, and 3 in fig. 6). 
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i) K = 0.086 de^deg/secj 
(Do = 13*65 radians per 
second; ^ = 0.057^. 


(b) K = 0.086 deg/deg/sec; 
(Do = 5*85 radians per 
second; ^ = 0.26. 


i<igure y.- Comparison of motions -with, two different autopilots idiose 
characteristics lie on a c\irve for Tpy'2 “= “ (corresponding to 

points 4 and 5 la fig* 6). 
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Figure 10.- Variation of optimum points in the ^cuo-plane for varying 

autopilot gearing ratio. 
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Figure l 4 .- Values of autopilot parameterB K and for vhlch an 

oscillation will occur in the lateral motion with damping slightly 
greater than that of the natiiral Dutch roll oecillation. C = 0 . 3 . 
Values of o) in ladians per second are given at representative points. 
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2 4 O Q /O /Z /4 /6 /Q 20 22 24 

Ol)^ , radians /6€C 

Figure 15,- Values of autopilot parameters K and 03^5 for which an 
oscillation will QCCVlir ^r\ 'h'hp T 1 -! 4 -V. m ( f\r\ 

— WJ.UJU -^1^2 " J-.uu btft;ona, 

the meilmum damping for K < 0 with ^ = O.J (the value of ^ in 
"this figure), Valu 6 B of co in radians par second ar© given ai 
representative points. 
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Figure 16 .- Values of autopilot parameters K and cOq for which an 

oscillation will occur In the lateral motion with ~ 0*60 second. 

^ =0.3. Values of 01 In radians per second are given at representative on 

points. ' 
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Figure 20 .- Comparison of the houndary of equal roots, the airplane- 
damping region, and a typical increased-damping region in the 
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Figure 20.- Continuei 
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(c) ^-0.9. 

Figure 20.- Concluded. 
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(a) Case I. 


(b) Case II. 


(c) Case III. 


Figure 23.- Calculated motions for three flight conditions of the airplane 
subsequent to a 5*^ sideslip disturbance. Autopilot characteristics: 

K = 0.l4 degy/deg/secj ^ = 0.523j 0\) = 9>49 radians per second. 
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(a) Case I. 


(b) Case II. 


(c) Case III. 


Figure 2k. - Calculated motionB for three flight conditions of the airplane 
subsequent to a 5° sideslip disturbance. Autopilot choracterlstlcB: 

K 0.12 degydeg/eecj ^ = 0.485j cDq = 8.8l radians per second. 
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(a) Case I. ("b) Case II. (c) Case III. 


Figure 25.- Calculated motlon6 for three flight conditions of the airplane 
subsequent to a 5° sideslip disturbance. Autopilot characteristics: 

K = 0.0T5 deg/deg/secj ^ = 0.389} oJc, = 7.40 radians per second. 
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